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We analyze Aharonov-Bohm radiation of charged fermions from oscillating solenoids and cosmic 
strings. We find that the angular pattern of the radiation has features that differ significantly from 
that for bosons. For example, fermionic radiation in the lowest harmonic is approximately isotropi- 
cally distributed around an oscillating solenoid, whereas for bosons the radiation is dipolar. We also 
investigate the spin polarization of the emitted fermion-antifermion pair. Fermionic radiation from 
kinks and cusps on cosmic strings is shown to depend linearly on the ultraviolet cut-off, suggesting 
strong emission at an energy scale comparable to the string energy scale. 
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I. INTRODUCTION 

The Aharonov-Bohm (AB) interaction [l|, Q between 
charged particles and thin magnetic fluxes is of much 
interest as it provides a physical consequence of a pure 
gauge field with vanishing field strength but non-trivial 
topology. Further, the physical effects emerge only in 
quantum theory and hence the AB interaction provides 
an example of a quantum, topological interaction. 

The classic Aharonov-Bohm setup involves scattering 
an electron off a solenoid, with non-trivial scattering ob- 
tained even for an arbitrarily thin solenoid, whereby the 
electron is exclusively localized in a region of vanish- 
ing magnetic field. A novel feature of the scattering 
is a periodic dependence of the scattering cross-section 
on the magnetic flux through the solenoid. If $ de- 
notes the magnetic flux in the solenoid and e the elec- 
tron charge, the cross-section is proportional to sin 2 (7re) 
where e = e$/27r. 

The classic AB setup was recently extended in an- 
other direction Q, where it was shown that an oscil- 
lating solenoid in vacuum can produce charged particle- 
antiparticle bosons from the vacuum due to the AB inter- 
action. The AB radiation rate also has the characteristic 
sin 2 (7re) dependence on the magnetic flux. 

AB radiation is relevant to the evolution of cosmic 
strings, which are similar to solenoids, except the mag- 
netic flux within them is massive, unlike electromagnetic 
fluxes in laboratory solenoids. Moreover, the gravita- 
tional analog of the AB effect can cause cosmic strings to 
emit light, even if the fields composing the cosmic string 
are unrelated to electromagnetic fields [1, H| . 

In this paper we will investigate fermionic AB radia- 
tion. One motivation is that the electron is a fermion. 
Hence fermionic AB radiation is what is relevant to os- 
cillating solenoids. The investigation is also relevant, for 
instance, to neutrino emission from cosmic strings by the 
AB process. A second motivation is that the spin of the 
fermion adds another degree of freedom to the emission 
and the polarization properties of the radiation are of 
interest. 

Our results show a significant difference between AB 
radiation of bosons and fermions. For example, if a 



solenoid aligned with the z-axis oscillates along the x- 
direction, the angular distribution of bosonic AB radia- 
tion is peaked in the y-direction. Fermionic AB radiation, 
however, is (approximately) isotropically distributed. 

The outline of the paper is as follows. In Sees. [TT1 
through IIV1 we use conventional interaction picture per- 
turbation theory to calculate the fermion-antifermion 
pair production rate in the small AB phase (e) limit. We 
consider AB radiation from an infinite straight solenoid 
oscillating perpendicular to its length in Sec. IIII1 a cos- 
mic string loop with kinks in Sec. IIV A[ and a cosmic 
string loop with cusps in Sec. IIVB1 In Sec. [V] we solve 
the problem using a different technique that does not as- 
sume that e is small. We call this the "moving frames" 
scheme and use it to obtain the sin(7re) dependence of 
the radiation on the AB phase, provided the motion of 
the solenoid is slow. We conclude in Scc. lVII and describe 
our conventions in Appendix [Al 



II. AB FERMION PAIR PRODUCTION 
A. Setup 

The interaction of fermions with the gauge potential of 
the thin solenoid or string is contained within the Dirac 
action 



= d 4 x ip {ijj) — to) ip. 



(1) 



(See Appendix lAl for conventions.) The relevant interac- 
tion term is 



(2) 



where is the classical solution around the flux tube 



2] 



with 



$ la 



(3) 



S a0 (x) = J drda^e ab d a X a d b X^5^ (x - X(a,r)) 
= J drda{X a X Si ' - X^X a ')S {4 \x - X(a, r)) (4) 
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FIG. 1: Feynman diagram for fermion-antifermion (which we 
have named e + e~ for convenience) pair production from a 
magnetic flux tube in motion. The black dot represents the 
classical gauge field of the flux tube and the solid lines are 
the outgoing electron and positron states. 



and X 1 - 1 (a, r) gives the position of the flux tube in terms 
of world-sheet coordinates a and r. 

We will need the Fourier transform of and this is 



given by 



,<I> 

■ 1 

2 



- -i— f a—^P 



(5) 



with 



S afi {k) = / d A x e +lk - x S af3 {x) 



dt / da 



~di da~' 



(G) 



where t — t and superscripts within square brackets are 
antisymmetrized . 

In the case of a straight solenoid, we will impose the 
dynamics by hand and consider oscillatory motion. In 
the case of a cosmic string, the dynamics will be given 
by the Nambu-Goto action as discussed in Sec. IIVI 



B. Pair production amplitude 

The amplitude for pair production is then given by the 
Feynman diagram in Fig. [TJ 



M(0 



~.S -,!/„, S 



J v {k,s;k',s')= U ±L V ]— , p^ = k fl + k^ 



where the s and k are the spin and momentum labels for 
the outgoing fcrmion state and the primed labels are for 
the anti-fermion. 

It will turn out, for all three cases of interest in this 
paper, S a @ can be factorized into an antisymmetrized 
product of two independent integrals, 



Furthermore, as can be checked explicitly, the I± and 
electromagnetic current are conserved 1 

P M i£=P lt J ,t = o 

This allows us to re-write 

e l *u a prPJ v S c # = P -3 ■ (1+ x I_) 
Po 



so that the amplitude now reads 

x I_ 



A4 = ^J-(I H 

Po 



e<I> 
2^ 



(7) 



Total rate: Our normalization of the Dirac spinors 
[3|) is such that the square of the amplitude itself is 
the total number of fermion-antifermion pairs produced. 
First let us evaluate the total rate, and hence sum over 
the spins of the final particles. In squaring the amplitude, 
we may then exploit the spin sums 



v%>v s k , -m, 



and the Clifford algebra {7^,7"} = 2?7 MI/ to deduce, 
k^k v + k' v k» ~ (m 2 + k ■ k')rf v 



knk 



K Q 



Within the small AB phase approximation, the result of 
this spin sum accounts for the entire difference in an- 
gular distribution for the outgoing fermion-antifermion 
pairs from that of the scalar case in 0. We see that, 
even at the level of unpolarized rates, the spin of the 
particles interacting with the vector potential give rise 
to significant observational signatures. 

For comparison, in the bosonic case the corresponding 
quantity is [3[ 



{k-k'Y{k-k') 1 

k k' a 



The difference in this expression versus the expression for 
the fermionic case gives rise to different angular distribu- 
tions for the AB radiation of bosons and fermions. In 
particular, AB bosonic radiation vanishes when k^ = k'^ 
whereas fermionic radiation does not. 

Now the square of the amplitude, summed over the 
possible spins of the outgoing particles, reads 



E 

s.s 



>' 2 =1S{ II+X1 - I2(m2+ '-' 



((1+ xl_)-k (1+ x I_)* ■k' + c.c. 



(8) 



1 For /J, see Eq. (|33}, and more explicitly, (TBJ, J37} and gOj. The 
identity for follows from the free massive Dirac equations, 
+ m)vf', = = ul (# - m). 
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with the "c.c." representing complex conjugation of the 
term preceding it. 

The x I_| will contain an infinite series involving 
the square of (5-functions, of the form (5(po — £S7)) 2 where 
f2 is the characteristic frequency of oscillation and £ is the 
radiation harmonic. This may be interpreted as a single 
(^-function multiplied by an infinite constant correspond- 
ing to the total duration of time divided by 2tt: 



(5(p -m)) 2 ->s(po-m)£- 



(9) 



where T is the total time duration over which the radia- 
tion is calculated. This can be justified formally by using 
the integral representation of the <5-function and setting 
the exponential to unity, 



6(0) 



T— ¥00 a— >0 J _t 12 27T 



T/ 2 



dt 



lim lim / — e lat = — 



T 
2^ 



The rate of pair production is then given by 
J2 S s ' l-^| 2 /^7 integrated over all kinematically possible 
momenta of the outgoing particles 



dN 
~dt 



(2tt)V 
T 



d 3 k 
(2tt) 3 



d 3 k' 



1 



lU x I_| 2 (m 2 



(2tt)3 plk k' Q 
k ■ k') 



(10) 



((1+ xL)-k (1+ x I_)* -k' + c.c. 



Spin- dependence: If one wishes to evaluate the spin- 
dependence of the scattering it is convenient to work 
out the amplitudes for the creation of particles and anti- 
particles with definite helicity. In other words the parti- 
cle of momentum k is assumed to have its spin aligned or 
anti-aligned with the momentum k corresponding to pos- 
itive or negative helicity. This is accomplished by evalu- 
ating the current J(k, s; k', s') taking the spinors u s k and 
vt, to be of definite helicity (see Appendix IX)) . An ex- 
ample of such a spin-dependent pair production rate is 
given in Sec. IIIII 

The spin-dependence of pair production has a general 
feature that can be deduced immediately from the form of 
the matrix element M. . Let us suppose that 1+ x I_ oc rh 
where rh is a unit vector. This assumption will be seen 
to be valid below for the cases of a straight oscillating 
string and for degenerate kinky loops discussed below in 
Sees. Mil and IIV Al For this case one can show that 



.VI xJm CSV,:i,- ' £S E -<7& 



(11) 



Here 



c c 

Ss,o = 2im ■ (n x h) sinh — sinh — , 

c c 

Sy; = 2 cosh - cosh — rh 
2 2 



(12) 



C C r 

+2 sinh — sinh — 
2 2 



(h • rh)n' + (n' • rh)h — (h ■ n')rh 



The rapidity is defined via cosh£ = k°/m, cosh£' = 
k'°/m and h and n' are unit vectors along the directions 
of the momenta k and k' respectively. Eq. (|11[) is a simple 
consequence of the forms of the spinors u s k and Vy given 
in Appendix [X] Here £ s is a spinor that corresponds to 
the spin of the particle in its rest frame; £ s < to the spin 
of the anti-particle in its rest frame. Thus we see that if 
the particle is measured to be up along the direction Ss 
in its rest frame, the anti-particle will definitely be down 
along the same direction in its rest frame and vice versa. 
To see this explicitly, if (J = (1,0) then the amplitude 
is maximized by taking Q = (1,0). However this choice 
of spinors means that the particle is spin up and the 
antiparticle is spin down. Hence there is a definite anti- 
correlation in the spin of the produced particles. 



III. INFINITE, STRAIGHT SOLENOID 

In this section we will consider the pair production due 
to an infinite, straight solenoid aligned parallel to the 
z-axis, moving in a sinusoidal fashion along the x-axis. 
Hence, with a = z, 



X a (t,z) = (i,£(t),0,z) 

= ^sin(Oi), Q > 



(13) 



Putting (fTBf into (|6|) we find that the t- and z- integrals 
may be factorized (into /+ and I- respectively) . To eval- 
uate 1+ we express the cos(f2f) in its integrand in terms of 
exponentials, and perform a cylindrical wave expansion 
via 



+00 



Jjp cos 6 



(14) 



The resulting Besscl functions in I + can be combined 
using the recursion relation 



V -JM = -(J v -i(z) + J v+1 (z)). 



(15) 



The integral for i_ simply gives a ^-function. Together, 



,Po 

Px 



(16) 



e=-co ^ 
it = S^2nS( Pz ) 

and therefore 

1+ x I_ = y(2^) 2 

OO 

\ i £ / p x 

where we have used 

J_ n (z) = (-l) n J n (z), n& {0,±1,±2,...} (17) 
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Making use of the formal identity 



L 

2tt' 



where L is the total length of the solenoid we find that the 
fermion-antifcrmion pair production rate (|10[) per unit 
length of the solenoid in sinusoidal motion is 



n> = jr / d 3 k J d 3 k r 

C—l s,s'—z 



where the differential rate 



d 7 N(ks,k's') 
dtd 3 kd 3 k' 



(18) 



0(vq £+2 ). (For I ^> 1, the Je(z) becomes exponentially 
suppressed; see equation ([2"5]).) Therefore, provided that 
SI > 2to, the first harmonic is the dominant production 
channel in the non-relativistic limit. 

We first, however, begin by performing a consistency 
check of our calculations based on the £ = 1 term. Using 
the pair production rate in the form (|10[) . and keeping 
only the leading order term in the series expansion of J\ 
and eliminating the k x k' x and k y k' y term by integrating 
over the appropriate angular coordinates, 



Af(ks,k's') 2 



(p x ^) s( Pz )S(p°-m) 



d 7 N(ks,k's') _ 
dtd 3 kd 3 k' ~ (2ir) 2 4k k' Q p 2 x 

(19) 

and s and s' represent the helicity of the particle and 
anti-particle respectively. The quantity TV is given by 



Af(ks,k's') 



(kok' + to 2 — ss'|k||k'| 



SS ill 

|k||k'| ^ z ^ z ^ x ^ x ~ kyhy. 



(20) 



Note that in the limit that the particle and anti-particle 
are ultra-relativistic Af(ks,k's') is negligible for s = s'; 
thus pairs are predominantly produced with opposite he- 
licity in this limit. The pair production rate summed over 
final state polarization is given by the simpler expression 



Af(k,k') = J2 N(ks,k's') 



(21) 



s,s' — ± 



= 4(m 2 + k k' - k z k' z - k x k' x + kyk' y ) 
^2([mf-{k x + k' x f -(k y -k' y ) 2 ) 

The second line in (f2"Tj) is the direct consequence of the 
term in curly brackets in (|10p . The third line, which 
is ^-independent, has been obtained from the first by 
imposing k z = —k' z and re-writing kok' = (1/2)(£S7) 2 — 
(1/2)/sq — (l/2)fcg 2 , using the constraints implied by the 
(^-functions. 

Note that when we square the amplitude we obtain a 
double sum over harmonics. This collapses to a single 
sum because 6(po — £Q)5(po — t'£l) is zero unless t = £' , 
since the requirement that both the (5-function arguments 
be null cannot otherwise be satisfied. Since the sum of 
the positive energies po = ko + k' cannot be zero or neg- 
ative, we have also removed all the £ < terms in the 
summation in Eq. (|18l) . 

Non-relativistic limit: In the non-relativistic limit, 
vq <C 1, we may utilize 



J ^-rFTT)©^ 1 + ^ 2 ))' I 



z|«l (22) 



to see that the contribution to the pair production rate 
at each harmonic begins at 0{vq 1 ) plus corrections of 



Eq HID 



dkk / dk'k' 
Ja 



dk 2 



x 5(k z + k' z )S{k + k' -Q) 
x (to 2 + k k' + k 2 z ) 



dk' 



f 2 v 2 

k k' in 2 



(23) 



This will match (|63[) from the small flux limit of the mov- 
ing frames perturbation theory calculation below. 

Now, the differential pair production rate is the inte- 
grand in (fTg)) . Upon integrating k z (or k' z ), the third 
line in (f2"Tj) indicates that the only k z = —k' z dependence 
in the integrand occurs in the kok' . Enforcing the con- 
straint ko + k' = by introducing a Lagrange multiplier 
we see that the integrand, and hence the emission rate, 
is maximum at k z — 0. (This is true of the full rel- 
ativistic emission.) Applying (|22|) to I = 1 and going 

, k y ) = k±(cos9,sin9) and 

of (ED then yields 



to cylindrical coordinates, (k x ,k y ) 

(k' x ,k' y ) = k' ± (cos 8', sin 9'), together with the second line 



d 5 N 



dkidk' d9d9'dk z 



= l,k z = k' z =Q) 



x (to 2 + k k' - k ± k' ± cos(6» + 9')) 



with 



ko 



k 2 



ft n 



The differential rate is largest in the xy-plane and when 
the outgoing particles' azimuthal angles are supplemen- 
tary, 9 + 9' = n. In contrast, boson emission is maximum 
along 2 9 + 9' = 2tt and dipolar (see Fig. [2]). 

We may proceed to employ the second line of (f2"Tj) to 
obtain the total rate of production of pairs with energy 
51 per unit length of the infinite flux tube 



N'{£ = 1) = 



e (Vo_y 
2(2tt) 2 \2flJ 



2 This corrects where it was mistakenly stated that the maxi- 
mum emission is along 8 + 8' = it. 




FIG. 2: The radiated power as a function of 6 and 0' . The top 
figure is for fermion radiation 1=1 with |fc| = IQ/2 = |fc'|, 
and the lower figure is for bosonic radiation with the same 
parameters. There is maximum radiation along the line 
9 + 9' = 7r in the fermionic case and the radiation is (approxi- 
mately) circularly symmetric. The bosonic emission is dipolar 
and along the 9 + 9' = 2-k line (plot taken from Ref. ||). 



£1 — m 



dkr 



dk z I dO d& 

Jo Jo 



[n 2 - (k x + k' x f - (k y - k') 



where 



(kx i ky) 
(k x i ky) 



k 



i 2 — k 2 (cos 8, sin( 



= y/(Sl - k Q ) 2 - m 2 - k 2 (cos 9', sin 9') 



In terms of the mass-to-energy ratio r = 2m j 'CI, the result 
is 



N'(£ 



v 2 e 2 n 2 

384 



6r 4 In 



2X/T - 7[16 - t{t(3t + 2) + 8}] 
+ 0(vt) 



(24) 



This pair production rate begins at v^e 2 Cl 2 /12 for r = 
and, for increasing t, decreases monotonically to zero at 
the threshold r = 1. As the fermion- antifermion pair 
gets heavier, they become harder to produce. 



Mode number ( 

FIG. 3: The radiated power as a function of harmonic for 
vq = 0.001 (squares), 0.1 (circles) and 1 (triangles). 



At very high harmonics (£ ^ 1), we may re-scale all 
momenta by £CI/vq and denote the new variables by over- 
bars, e.g. k = kvo/(£Q). Then we invoke @ 



cosh< 



exp (— v(a — tanha)) 
\j2itv tanhev 



a > 0, v > 1, 
(25) 



to say that, within the integration limits < p x < vq <C 
1, 

/ Je(£px)\ 2 exp[— 2£(ao — tanhao) + 2 ln(cosha )] 
\ Px J V27r£tanh ao 

(26) 

where ao is the positive solution to the equation 
coshcio = 

The exponent on the right hand side of the asymptotic 
expression (|26[) is a monotonically increasing function of 
px, and takes on large negative values for small p x (i.e. 
large ao), while tanhcto ~ 1. Hence the pair production 
rate (fT5|) is exponentially small for £ ^> 1, in the non- 
rclativistic limit. In Fig. [3] we show the power emitted 
per harmonic for a few values of vq. 

High frequency limit: When CI 2m, such that 
the fermion-antifermion pairs are produced with very 
large momentum, to a good approximation, they may 
be treated as effectively massless. If we re-scale all mo- 
menta by the energy, k = £Clk etc., and integrate over 
k',, the emission rate in the xy-plane is 



d^N 



dk ± d0d0'dk 



■(k z = k' z = 0,m = 0) 



£Cle Je(£v px) 



2tt 



k±(l-k±)(l 



<e + e')) 



with 



k± cos 9 + (1 — fcjj cos9' 



The rescaled momentum, k±, lies in (0,1) and so p x 
lies in the interval (cos 9, cos 9') for cos 9 < cos 9', and 
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(cos 8', cos 8) for cos# > cos#'. Therefore the absolute 
value of the argument of Jg is always less than its order, 
I^oPkI < £ v o < £■ For I = 1, as a first approximation, we 
may recall (f2"2"j) , and observe that (Ji(£vop x )/p x ) 2 stays 
roughly constant within the range of interest at hand. 
Hence, the maximum emission occurs for kj_ m 1/2 and 
8+8' = tt. For I > 2, we may employ (jT5j) followed by the 
fact that both the first turning point and zero of Je(z) 
occurs only for z > £, to argue that (Je(£vop x )/p x ) 2 is a 
monotonically increasing function of p x — this means the 
kj_ at which there is maximum emission now shifts away 
from 1/2, with the direction depending on whether cos 9 
is greater or less than cos 9'. 

At higher harmonics, I 3> 1, Ji becomes exponentially 
suppressed, according to (|25|) and (|26]l. in most of the 
{6,6') plane. If the motion is relativistic (vq ~ 1), how- 
ever, the asymptotic formula [H 



Nambu-Goto action 



Jt(z) 



1 /sin(7r/3)r(l/3) 
Sir 



(z/6)V3 

sin(27r/3)r(2/3) 
(z/6) 2 / 3 



(z - + 



(27) 



»1, 



valid for z ~ £, says that (Jg(£vop x )/p x ) 2 will transition 
to an inverse power law in £ and hence peak in the region 
where \p x \ lies closest to unity, namely, where cos 6? = 
cos0' = ±1. (Note that ( J e (-z)/(-z)) 2 = (J e (z)/z) 2 .) 
Therefore, ignoring the (1— cos(9+8')) factor for now, the 
peak occurs at (9, 9') = (0, 0) or (tt, tt). The shape of the 
peak is determined by the constant p x contour lines on 
the (6,0') plane near (0,0) and (it, it). They are ellipses 
with fc^-dependent eccentricity, because 

p x (0, 0) = 1 - \(k^9 2 + (1 - k ± )6' 2 ) + ... 



Px(k,1t) = -1 + -(fcj_( 



- tt) 2 + (1 - k x )(8' - tt) 2 ) + 



Now we include the effect of the (1 — cos((9 + 8')) factor 
which vanishes along the lines 6 + 6' = lirn, where neZ. 
The peaks at (0, 0) and (tt, tt) lie precisely on these lines. 
So the multiplicative factor (1 — cos(8 + 8')) modulates 
the peak and splits it into two, one on either side of the 
6 + 8' = 27m lines. 

To summarize, at large oscillation frequencies and rel- 
ativistic solenoid speeds, we expect the high energy spec- 
trum of the fermion-antifermion pairs to be produced pre- 
dominantly in the xy-plane, in narrow pairs of nearly 
collinear beams, in the ±x directions. This beaming 
effect is due to the presence of the Bessel functions 
Ji(£vop x ), common to both the fermion and boson am- 
plitudes, and hence will occur in the boson case too Q. 

IV. COSMIC STRING LOOPS 

We now turn to fermion pair production from cos- 
mic string loops for which the dynamics is given by the 



Sng = — A 1 dr da 



(28) 



where the world-sheet metric is 

lab = V^daX^d b X" , 7 = det( 7ab ) (29) 

and £ a = (r, a) are world-sheet coordinates. 

Referring to the Nambu-Goto dynamics for relativistic 
strings in (|28p , we will choose a such that constant t- and 
cr-lincs on the string world sheet are orthogonal, and the 
energy per unit (proper) length of the string is constant, 



X' ■ X = 0, 



X' 2 + X 2 = 0. 



(30) 



These choice of coordinate conditions (|30|) put the world 
sheet metric 7 j, into a conformally flat form, and lead us 
from the Nambu-Goto action ([28]) to the wave equation 
for X^ 



(d?~d 2 )X»(t,a)=0 



(31) 



The X^ solution may be written as an average of lcft- 
and right-movers, L M = L»(o + t) and = i?^(cr - t) 
respectively, 

where we also have introduced world sheet light cone co- 
ordinates 



a± = a ± t, 



and 



L° = a 4 



The world sheet coordinate conditions (|3T)]) translate into 
constraints on the spatial components of X^, 



R' ; 



1 



We will focus our attention on two specific loop con- 
figurations. The kinky loop has tangent vectors to the 
string that are discontinuous at isolated points; whereas 
the cuspy loop has isolated points that reach the speed of 
light periodically. These configurations satisfy the equa- 
tions of motion (|3"Tj) that follow from ([2"5]l . 

If the coordinate length of a given cosmic string loop 
is L, the a- integration limits for its corresponding S^, 
(see Eq. ([6])) runs from to L. Since we have closed 
loops, X and derivatives dtX^ and daX^ are periodic in 
a as well as t. To facilitate the computation, we utilize a 
formal device introduced in Q that allows us to factorize 
into a product of two one-dimensional integrals, one 
for each world sheet light cone coordinate. This is based 
on the observation that, for a periodic function / with 
period L, the integral over one period of / is equivalent 
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to the integral over the real line, divided by the (infinite) 
number of times the former has been over-counted 



i r+oo 

o dz(0) 



where 



5 Z (0) = ^ — , £eZ 



la/L 



We may now extend the a integral in S^ v to the entire 
real line using (|32p . before changing variables from (t, a) 
to (a + , cr_ ). S^ v then factorizes into 



guv _ _j\mj"] 

r> + 00 



(33) 



1 

I+ = - dcr + d + L a e ip - L/2 

2 J-oo 

7 - = ^T7^ / da-d-ITe**-*' 2 , 
2<5 Z (0) J_ OQ 



where the derivatives are with respect to a± 3 . Periodicity 
in a and t implies that the integrands in (|33j) are periodic 
in a + and er_ and we can replace them in Eq. ^ with 
their discrete Fourier series expansions, 



l=-oo ^ 

x ( L ^±d+L a e- ll2 ™+/ L e-*P 
Jo L 



L/2 



(34) 



r = 



2tt 

W) 



da- 



Jo L 



The 1+ x I that follows from (|34|) will be an infinite 
sum involving (5(po — £4w/L)) 2 . We write 

(5(p - £4ir/L)) 2 = S( Pa - £4n/L)d z {0)-^ (35) 

47T 

and the 5 Z (0) in (J35]) will cancel that in J_ (|34|) . 



3 If we had not applied 11321 , upon converting to light cone coor- 
dinates, the inner <r±-integration would have limits that depend 
on the outer integration variable, and S^ 1 " would not factorize. 



A. Kinky Loops 

The "degenerate" kinky loop solution we will consider 



R(o-_) 



(7 i A 



< a + < f 



(L - cr+)A I < CT+ < L 



B 



< CT- < f 



(X - cr_)B 4 < <r_ < L 



(36) 



where A and B are unit vectors. This loop is degenerate 
because it consists of four straight segments and is kinky 
because of its four corners. The four straight segments 
propagate with constant speed but shrink and expand 
due to the motion of the kinks. 
Now, denoting 

p A =p t A\ p B =PiB\ 



and putting the kinky loop trajectory (|36f into (|34[) then 
leads to 



n= E (PA,-p A)s( P0 -^£ 

( 



167re i(PAi/8-^/2) n 

T, j 

2 



1 00 / A \ 

x lfafaW '4^ + ^] (37) 



and 



16 E 



,4^ gi(PA+PB)i/8 



x sin — L 1 sin 

V 8 2 / V 

xs(p - A x B 

The power radiated from the kinky loop due to fermion 
pairs emitted with energy £fl, for a fixed £, is the corre- 
sponding pair production rate multiplied by the energy 
4tt£/L, 

^(K)_ 9 2 / 647r\ 2 /47rA 3 f d 3 k d*k> 1 
1 ~ \ L ) \ L J J (2tt)3 ( 27 t)3 
sin 2 (paX/8 - ttI/2) sin 2 (p B L/8 + ttI/2) 

X (pI-pIWb-pD 2 

x [(A x B) 2 (m 2 + k ■ k') + 2A x B • k A x B • k'] 
x 5 (p - *Et) (38) 



The total power radiated is then 

N 
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FIG. 4: Plot of power radiated in the form of massless 
fcrmion-antifermion pairs in units of 64e 2 /(-7r 3 L 2 ), from the 
degenerate kinky loop with A and B perpendicular, as a func- 
tion of the cut-off mode number N. 



We have truncated the summation at some large inte- 
ger N because it will turn out, just like in the bosonic 
case 0], that the total power obtained from summing to 
£ = oo will diverge. The cut-off is related to the round- 
ing off of the kink and may be estimated as the ratio 
of the length of the cosmic string loop, L, to its width 
w. Taking L ~ 1 Mpc and w ~ 1 TeV -1 , we obtain 
N ~ 10 41 . For the electron - the lightest electrically 
charged fermion - note that the relevant range of mode 
number is 10 35 < £ e + e - ^ 10 41 for such an L and w, 
where the lower limit is determined by the product of 
the electron mass and L. In the high energy range of the 
spectrum, which gives the dominant contribution to the 
power emitted, the fermions can be treated as effectively 
massless. 

When m = 0, by re-scaling the momenta in the integral 
([55)1 via (k, k') = (4w£/L)(k, k'), the sole dependence on 
the summation index I occurs in the trignometric func- 
tions, which in turn can be summed using the formula 



N 



^sin 2 (» sm 2 (£y) 



M sin(Afx) sin(My) 



> sm x 



i sin y 



sm{M{x-y)) + sin(Af (x + y)) M = m + 1 
16 sin(.T — y) 16 sin(x + y) ' 



The resulting integral was then evaluated numerically 
with Mathematica. 

In Fig. [4] we have plotted the total power E(N) as a 
function of the cut-off N, for m = and for the "square" 
loop which has A • B = 0, from N = 1 up to N = 10 50 . 
We see that the power radiated per mode £ is indeed 
independent of £, in this massless limit, and hence the 
total power emitted grows linearly with N ~ L/w and 
can be very large. 



B. Cuspy Loops 



The cuspy loop solution we will consider is 
L(a + ) = - sm — j,0,-cos f_ 



™, \ L ( . (2-kv- 



1 ~L~ ' 



The loop has cusps because there are points such that 
<9+L = — cLR. These occur at 2ira±/L = 0,7r at which 
point the velocity of the string is ±x, i.e. the point on 
the string reaches the speed of light. 

Putting the cuspy loop trajectory (|39|) into (|34[) . com- 
bining the sines and cosines occurring in the exponen- 
tial into one trigonometric function before performing a 
cylindrical wave expansion with (| 14|) then gives us 



2tt 



— — OO 

oo 



(40) 



with 



<P(+\e)(Px,Pz) = i l Ji (J^VpI +Pl S j CX P ^arctan (jf^j 
f(-\i){Px,Py) = -i e Ji C^\/Pl+Pyj CX P ( il 



: arctan 



Odd 
dp x ' dp y ' dp z 



as well as 



j j 16tt 3 ^ / 4tt 
1+ x I_ = — 2^ s[po-—e)nt 



with 



n; = y x 9 P ^ (+ |_£) j x (z x <9 p </?(_ 



The derivatives on Jg occuring in 1± and may be car- 
ried out with the aid of 

d,J v {z) = ^(J v -i(z)-J v+ i(z)) (41) 

The power radiated from the cuspy loop due to fermion 
pairs emitted with energy £Q, for a fixed £, is the corre- 
sponding pair production rate multiplied by the energy 
4w£/L, 



E 



(C) 



d 6 k I d 3 k' 



2 e 2 



£Lk kl 5 [ P0 



4n 



x <{ \n e \ 2 (m 2 + k- k') + (nt • k n* e • k' + c.c.) } (42) 
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FIG. 5: Plot of power radiated in the form of massless 
fcrmion-antifermion pairs in units of e 2 /(2nL 2 ), from the 
cuspy loop, as a function of the mode number I. 



Using Mathematica, we evaluated (|4"!2)) numerically for 
I = 1 through I = 1000 and m = 0. Fig. [5] displays the 
resulting power radiated in the form of massless fermions 
at each harmonic, with energy Att£/L. This provides evi- 
dence that, at large mode numbers (I ^> 1), the power ra- 
diated varies very slowly with £; though a thorough anal- 
ysis would have to employ more sophisticated numerical 
techniques (or semi-analytic ones, using (|25)) and ([27)1 ) 
to evaluate (|42j) for the astrophysically relevant range of 
10 30 < i < 10 50 . 



V. INFINITE, STRAIGHT SOLENOID VIA 
MOVING FRAMES PERTURBATION THEORY 



In the second half of this paper, we provide an alternate 
calculation of the fcrmion pair production rate from the 
infinite straight solenoid, introduced in Sec. IIII1 aligned 
with the z-axis, moving in sinusoidal motion along the 
x-axis. The spatial location of the solenoid func- 
tion of time t is given by X(i, z) = (£(i),0,z), with 
£(<) = (vo/il) sin(f2i). This non-relativistic computation 
will only capture the first harmonic of the infinite sum 
obtained in Sec. lIIII - see (|44p below - but will make man- 
ifest the periodic dependence of the pair production rate 
on the AB phase e = e&/2ir that is expected from such a 
topological interaction. It will also serve as a consistency 
check on the results in Sec. Mil 

Let X(s|+) (x)e _lBs * denote the positive and 
X( s |-)(x)e +lBs * the negative energy solutions to 
the Dirac equation in the presence of the gauge potential 
of a static, infinite straight solenoid aligned along 
the z-axis centered at (x,y) = (0,0). (The subscripts s 
and s' in this section denote, collectively, all the rest of 
labels that come with the solutions.) We will expand the 
Dirac operator, ip, in terms of "shifted" mode functions, 
X(s| ±)( x ~ £i z)e TtEst . Then we can show, as derived in 
Sec. IV Al that the fermion-antifcrmion pair production 



amplitude for a solenoid moving along the x-axis is 
provided by the product of the integrals 

/+oo 
dt'at')e^' +E ^' 
-OG 

x J dV X { s|+) (x')^ XM -)(x'), (43) 

where the left-hand side of the equation is an inner prod- 
uct in the Heisenberg picture. 

If we specialize to the sinusoidal trajectory considered 
in Sec. IIII1 the time integral in (|4"3")l can be evaluated 
immediately to yield conservation of energy 



dt'£(t'y( E °' +E ° )t ' = -ittv 6(E s + E s , - Q) (44) 



(The 6(E S + E s i + ft) term was dropped because, for 
ft > 0, the sum of the two positive energies E s and E s i 
will never be negative.) The presence of the overall factor 
vq in (|4"4"| and the absence of any dependence on the 
trajectory in the volume integral in (|43"]) tells us that the 
perturbative scheme in this section is non-relativistic. 

We will proceed to solve for the complete set of modes, 
X(s|±)( x )e =F * £ ' st , evaluating (|43|) for necessary combina- 
tions of s and s' , before summing the squares of the re- 
sulting amplitudes to obtain the pair production rates 
(1521) and (Pl> . 



A. Derivation of Equation (|43|) 

We expand the Dirac operator, within the Heisenberg 
picture, in terms of X( s i±) (x — £, y, z)e TlEst 

s 

+ (3t(t) X( s\-)(x-^y,z)e + ^ t ] 

The a s (t) and /?J(i) are time-dependent operators that 
obey appropriate anti-commutation relations and ensure 
positivity of the Hamiltonian. They also satisfy the equa- 
tions of motion obtained from the Heisenberg equations 
of motion: 

idfip = hx>il> (45) 
/i D = eA Q - i^-fWj + 7717° 

For a solenoid at rest at £ = 0, the operators 

a s = a s (t ->• -oo), b s = fi s {t -co) 

acting on the zero particle state at t — > —oo destroy, 
respectively, a fermion (a s ) and an antifcrmion (b s ) asso- 
ciated with the wavefunctions X(s|±)( x )- I n particular, 

{a s ,al,} = 6 S>S >, {b s ,bl,} = 6 S>S ' 
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Then, to leading order in the Hamiltonian, the zero par- 
ticle to two particle amplitude is 



—i(0:t = — oo 



a s b S ' / dtd 3 x ip^idtip 



0; t = — oo 



(46) 



In (J46J) . we have re-cxprcsscd the Dirac Hamiltonian den- 
sity 

using its Heisenberg equations of motion (|4"5)l . (This is le- 
gitimate because the ip's appearing in the amplitude (|46[) 
are operator solutions to the Dirac equation.) The rea- 
son for replacing with id t is the following. Since the 
time evolution of the a s (t) and /3 s (t) away from a s and 
b s are determined by the interaction between ip and A„, 
to lowest order in interaction, we may now approximate 



a s {t) 



so that 



{a s ,^} = X( s \+)(x - ty,z)e tEA 
{b s , ip} = X( s |-) {x - f , y, z)e +lEst 

By anticommuting the a s and b s > in (|46|) to the right, 
and noting that - because the solenoid is moving solely 
in the x-direction - 

9tx(x - £(*), y, z) = -£(i)d x x(x - £(t),y, z), 

we then arrive at (|43[) after shifting integration variables 

x — > x — £. 



B. Modes with k z = 

The solution for the gauge field of an infinite straight 
solenoid along the z-axis with a = z has only two non- 
zero components 



F 0y = 6(x-S(t))6(y), 

F xy = -B z = -$ S(x - £(*))<%), 



(47) 



The gauge potential that yields (|4T|) is, in Cartesian co- 
ordinates, 



A„ 



(0,0,-$e(x-ti(t))6(y),0) 



(48) 



We wish to solve for the fermionic modes, x, m the sta- 
tionary solenoid (£(t) = 0) case. It is easier to do so using 
the Lorenz gauge, where in cylindrical coordinates, 



.4 



(S) 



|(S) 



4 s ) = 4 s ) = o. 



(49) 



where the superscript S denotes "stationary" . 



With the gauge potential determined, the Dirac equa- 
tion now reads 4 



■<eAg») 



-iEt 







(50) 



Translational symmetry along z implies that x has the 
form 



X(x) = g(r, 6)t 



(51) 



Following Alford and Wilczek @], to exploit the z 
translational symmetry, it helps to find a set of {7^} ma- 
trices such that all of them except j 3 are block diagonal. 
This way, in the reference frame where there is no mo- 
mentum along z, i.e. k z = 0, the Dirac equation splits 
into a pair of coupled equations, each involving only 2 
component spinors. The set of {7^} we will use here are 
defined relative to the ones in the chiral basis {7^} in 
(|AT|) as 7^ ee U^W with 



The new 7^ arc 





- 


—i 





i 




1 


-1 





-1 












-1 





-1 






i 





—2 







a 3 








'ia 1 





cr 3 




7 1 







-1a 1 



7 



—10 








7 



-il 
-il 



Now the k z = solutions are 

Xk(t,r,6,z) 



XA 

XB 



-iuj±_t 



EE ±y/tf 



(52) 



with 



(o-l) 



J tJl („ 1+ i +e )(fcr)e l9 

m+w ± ^i(«i+e)( fcr ) 

^J Mn2+1+e) (kr)e i6 

J a 2 {n 2 +e){kr) 



(53) 



where ai.2 = sgn(m,2 + e), ^1,2 € Z, e = e$/27r and the 
normalization factors A/a.b G C The spinors xa and xb 
are the two independent solutions for each set of positive 
or negative energy states and we have chosen our z-axis 
so that e = e$/(2?r) > 0. 

For 711,2 + e > and ni,2 + e < — 1) the signs o\ y % in the 
indices of the Bessel functions are fixed by the require- 
ment that the solutions Xk must be square normalizablc. 



4 In this section, the mass term in the Dirac equation has a + sign, 
as opposed to the — sign implied by JTJ . To convert solutions for 
one into solutions for the other, multiply the Dirac spinor by 75. 
Since this corresponds to a change of basis, this choice of sign for 
the mass term does not affect the results for the inner products. 
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Specifically, the radial integral must converge at its lower 
limit. It contains a factor of r, i.e. J Q drr, whereas the 
Bessel functions behave, for kr <C 1, as J v oc r v (see 

When n 12 If £ (— 1, 0), we have chosen <j 12 = —1 so 
that the solutions Xk remain square normalizable for all 
values of magnetic flux $. Even though both signs are 
allowed by square normalizability for $ =/= 0, our choice 
of sign gives solutions that, in the zero flux limit, join 
smoothly onto the cylindrical wave solutions to the non- 
interacting massive Dirac equation (i$ + m)ip = 0. 

It is worth expanding upon this ambiguity and the res- 
olution we have adopted here. This is a subtle issue dis- 
cussed first in context of strings in Ref. @ (see also Q , 
Q and 0). The basic problem is that for certain angu- 
lar momentum channels there are four solutions to the 
radial equation that are normalizable rather than two 
(for a given energy and z-momentum). Thus additional 
boundary conditions must be specified at the origin to 
determine the two physical solutions. The set of pos- 
sible boundary conditions is restricted by the condition 
that only those solutions are permissible that have zero 
radial current at the origin ("self-adjoint boundary con- 
ditions"). However this does not by itself fully specify the 
boundary conditions and some additional physical prin- 
ciple or regulation scheme must be invoked. In effect we 
have made a special choice of boundary conditions above. 
Our choice is natural and well-motivated for the follow- 
ing reasons: (1) It extrapolates smoothly to the case of 
integer flux. (2) The transition rate calculated using this 
boundary condition agrees with the perturbative result 
in their common domain of validity (see Sec. lVFl below). 
(3) One can imagine natural regulation schemes in which 
this boundary condition will arise Q- Still it must be 
kept in mind that the results would come out different if 
different boundary conditions were used. 

Shifting 2 in (|52p by an integer, ni t 2 — > rii^+m with 
m, G Z, takes us from one solution to another. Hence we 
can absorb the integer part of e in the label rii^ and, 
without loss of generality, choose the AB phase e to lie 
in the interval [0, 1). We will denote the fractional part 
of 6 as k. defined by 

k = e mod 1 

In terms of k, note that the set of mode functions here 
(and the set with k z ^ below) is now manifestly peri- 
odic in the AB phase e<f>. Hence the radiation rate that 
follow will also enjoy this periodicity. 



C. Modes with arbitrary k z 

Above we have found all modes with k z = 0. Now we 
perform boosts along the solenoid to obtain modes with 
k z 0. Denote the Lorentz boost by B and so 

Xk -> B ■ Xk 



Under the transformation 

t — > t cosh rj — z sinh r\ 
with the rapidity parameter rj defined through 



tanhr? = — 



(54) 



the boost matrix B is 



B=U-B c -U-\ 



B r = 



.77/2 1-<T 3 



-y/2 1+ct 3 







c n/2 l+a> + c -n/2 l-<r 



Denote xi as the boosted Xk with A/b = 0; and xn as 
the boosted Xk with A/a = 0. The two independent solu- 
tions for each set of positive and negative energy wave- 
functions, in the Lorenz gauge for (|4T)]) , are then 



X(i,n,k,k,\±)(t> r >9, z ) 

cosh(?7/2) J CT(fi+n+1) (kr)e lf> 
cosh(r)/2)aX±(k)J aiK+n) (kr) 
-i smh{j]/2)J a{K+n+1} (kr)e ie 
i sinh(?//2)crA± (k) J a ( K+n ) (kr) 

xN±(k,k z )e* ikot e in9 e ±ik * z , 
X(u,n\k>,k> z \±)(t,r,9,z) 

i S inh(r//2KA ± (k')Ja' ( K +n' + i) Wry 6 

-ismh(ri' /2)J (T ,( K+n , ) (k / r) 
cosh(7?'/2)a'A ± (k>)J a , {K+n , +1) (k'r)e ie 
cosh(?//2) J CT /( K+ „-) (k'r) 



(55) 



y.N±{k' ,k' z )e^ lk ' ot e m ' 



8 ±ik' z z 



now with 



1 if 7i, ri > 
-1 if ri, n' < -1 



and 



A±(fc) 



m + u>± 



, Af±(k,k z 



uj± + m 



2fcn 



(56) 



k, fc' > 0; k z ,k' z eR; n, ri = 0, ±1, ±2, . . . 
Wc have normalized our solutions using 

' dr r JJkr)JJk'r) = 6 ^t_^l 
Vkk 1 

k,k' > 0, Re(i/) > -1 

',ir 

dee i(n-n')e = 2TrS n , n ,, n,ri eZ 



+ oo 

dze z{q ~ q ' )z = 2nS{q - q'), q, q' e R 
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such that 



,fe,fc,k B )( x ')^(^'^'^'.^kE)( x ') 

2 S(k~k')S(k z - k' z ] 



jA,A'S n ,n'(2-!T)' 



Ikk' 



(57) 



(7 B ,(7^ =±, A, A' G {1,11} 



We will be comparing the results here to those in 
Sec. IIII1 where perturbation theory was performed with 
plane wave solutions written in Cartesian coordinates. 
To make the comparison, it is worthwhile to note that 
the normalization in ((57)) is consistent with the (k|k') = 
(27r) 3 <5 3 (k — k') one would otherwise have obtained if 
Cartesian coordinates were utilized. Moreover, 

(2^) 3 5 3 (k-k') 

= (27r) 3 (5(0 - (j)')S(k - k')S(k z - k' z ) (kk' f 1/2 . 

The \Jkk' in (|57[) is the (symmetrized) Jacobian when 
transforming from Cartesian to cylindrical coordinates. 
The (2ir)5(4> — 4>') is the completeness relation 



271-5(0- 0') 



+ 00 



e in (<f>-<f>'\ |0-0'|e [0,2tt), (58) 



Denoting our mode functions here by \n, k, k z ) (and sup- 
pressing the A and erg dependence), we see that ([58)) 
accounts for the missing (2ir) in ((571) because \n, k, k z ) is 
the nth term in the Fourier series expansion of |k), the 
mode functions written in Cartesian coordinates, since 

|k) = \k,k„<j>) =^|fc,fc z ,n)e™f 

n 

and hence (fS"71) yields 

Y^(k,k z ,n\k' ,k' z ,n')e m ' 4 '' e _in * = (k|k'). 



contain a e 4e(e mod 27r \ The e 4K < 9 mod 2?r) in ([59]) would 
then follow from the re-definition of the Fourier mode 
label n that introduced k. 

E. Matrix Elements for Moving Solenoid 



With the solutions ((59)) in hand, we are now ready to 
evaluate the integral (f4"3")) . It is only necessary to cal- 
culate (|43|) for positive and negative solutions from ((59)) 
with Fourier mode labels (n,n') = (—1,0) and (0,-1), 
where the n refers to the positive energy solution and n' 
to the negative energy solution. For other values of n and 
n', the factor of d x X( s >\-) in Eq. (j4"3")) yields the difference 
of two negative energy solutions in (f5"9")) , which by ((57)) 
has zero overlap with all its positive energy counterparts. 

For n' = 0, dxX^^ 1 ^ is a linear combination of the 
n' = — 1 and the n' = +1 negative energy solutions in 
(1591 . except the former has a 1 = +1. For n! = — 1, 



(<r'=-i) 
( s >\- 



is a linear combination of the n' 



-2 and 



nf = negative energy solutions in ((59")) . with the latter's 
cr' = — 1. In addition, the ml = — 1 solution with cr' = +1 
in c?a;X(n'=o| — ) remains orthogonal to all positive energy 
solutions except possibly for the n = — 1 case; likewise 
the n' = solution with a' = —1 in 9xX(n'=-i|-) is 
orthogonal to all X(s|+) except perhaps X(n=o|+)- 

As we will see in a moment, these two overlap integrals, 



d 3 a 



(<r=-i) 
X( n= _i|+) 



(<t'=+1) 

X(„'=_i|_) 



and 



d 3 x' 



(<7=+l) 

X(n=0|+) 



\ 



(<7'=-l) 

(n'=0|-) 



do yield non-zero answers and are therefore the only ones 
contributing to the pair production rate. 

Computing the volume integral in ()43)) for solutions 
with the Fourier modes given by (n, n') = (—1,0) and 
(0, —1) requires the integral 



D. Transformation to Axial Gauge 

The above modes for a stationary solenoid, £ = 0, have 
been found in Lorenz gauge. However, to find overlaps 
and the radiation rate, it is easier to work in axial gauge. 
Hence we now transform the modes to axial gauge with 
An in g8]), 



A M = (o,o,-$e(.T)%),o), 



(Cartesian) 



The x solutions differ from its Lorenz gauge counterpart 
in ((55)) only by a phase factor 



X(axial) 



Ak{9 mod 27rV, 

e X(Lorenz, Eq J55J) 



(59) 



In ((59")) , to see it is k and not e that should occur in the 
phase, refer to the solutions ((52")) before the introduction 
of k, and note that the axial gauge version of ()52l) would 



dr r J v (qr)J_ v (kr) 



Pr 



2 sin(7r^) 
■K{q 2 - k 2 ) 



(!)" 



C0S(7TZ/) 



S(k-q), 



k,q > 0, Re[u] > -1 (60) 



where Pr denotes the principal value. 
Define 



and 





The result then reads 
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d3x 'x\ A ^ Kkz \ +) {^)d x X(A>^M'K\-)^)) 



k^k' 



= ——5(k z + k' z ) sin(7TK)v / kk' 



' k' x K 



uj k + w k > 



Uk + Wfc' 



5^ A ,kj^^ + m,k' IUJk+m 



Lu k + m 



{uj k - uj k >)\/ujk + niy/ujk' - m 



uik> - m 



8 A 5 A ,kk' - 5 A 5 A ,(u>k + m)(uk' - m) 



5 A 6\,k\l^^ + 6WlkJ UJk '~ m 



tJk' - m 
H + -isgn(k' z ) 



(uj k - U3k')y/U3k + Vly/Uk' - m 



bj k + m 
i ill 



5 l A 8 A , (uj k + m)(uj k i — m) — S l A S l A ,kk' 



(61) 



Direct substitution of the Bessel integral Eq. ([60]) 
would lead to an additional term in the transition am- 
plitude above that is proportional to 5(k — k'). We have 
omitted this term following the discussion in Ref. Q. As 
explained there, the integrals must be carefully regulated 
in order to exclude spurious terms that lead to particle 
production even in the limit of zero flux. This is accom- 
plished by performing the integrals over a finite volume, 
imposing suitable boundary conditions, and then taking 
the infinite volume limit. Such an analysis eliminates the 
term proportional to S(k — k') (see Sec. IIIB of ref Q for 
more details). 

The expression in (|5T|) yields 8 different channels for 
pair production to occur. Moreover, at the level of indi- 
vidual amplitudes, the periodic dependence on the mag- 
netic flux is already manifest, since they depend on the 
AB phase only via k = e mod 1. 



F. Moving frames pair production rate 

Inserting (|61[) into (|43|) and summing the squares of the 
resulting individual amplitudes then gives us the rate of 
fcrmion-antifcrmion pair production per unit length of 
the infinite solenoid 



N' 



dk 7 



dk' 



dkk / dk'k' 
Jo 

xS(k + k' - n)S(k z + k' z ) 

X i£^M (m 2 + k l + k k' Q ) 



in^kok^ 



(62) 



In the limit when k <g. 1, the last term (in parenthesis) 
goes to 2 and sin 2 (7r«) w tt 2 k 2 , and 



,2 ,,.2 r°° 



N' 



An 2 



dkk / dk'k' / dk, / dk' 



x6{ko + k' -n)6(k z + k' z )[l + 



k 2 + m 2 s 



(63) 



This agrees with (|23|). the non-relativistic limit of the 
interaction picture perturbation theory result, if we iden- 
tify k e. 



VI. CONCLUSIONS 

We have solved for fermionic radiation first from os- 
cillating electromagnetic solenoids and then from cosmic 
string loops. For the solenoid we have done the calcula- 
tion in two different ways, first using a small AB phase 
approximation, and second by considering slowly moving 
solenoids. We have evaluated the angular distribution of 
the fermionic radiation from the solenoid, and the total 
power emitted from cosmic string loops and cusps. Our 
results can be compared to the results of Ref. [|| . 

The total power emitted in bosons and fermions is 
very comparable. For example, both are proportional to 
£ 2 Wq0 2 for the lowest harmonic of the oscillating solenoid 
(see Eq. ([24]) ). However, the angular distributions of the 
radiation in the two cases are quite distinct. To high- 
light the difference, we show the angular distribution in 
both cases for the lowest harmonic emission from an os- 
cillating solenoid in Fig. [5] We also find that the fermion 
and antifermion are preferably emitted in opposite helic- 
ity states and discuss the spin distribution. Just like in 
the bosonic case, fermionic AB radiation from kinks and 
cusps on cosmic strings is ultra-violet divergent for mass- 
less fermions, with a linear dependence on the cut-off. 
This may translate into a significant amount of radiation 
of neutrinos from strings with which neutrinos have an 
AB interaction. 

Our results ought to apply to the low energy end of the 
emission spectrum of electrons from idealized solenoids, 
where the wavelengths of the particle pairs are much 
longer than the diameter of the solenoid. A more re- 
alistic theoretical investigation would have to take into 
account the finite width of the solenoid itself. 
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In Ref. [3(, the gravitational analog of AB radiation 
was also discussed. Via the same analogy, we also ex- 
pect cosmic strings to radiate fcrmions. We leave that 
calculation for future work. 
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Appendix A: Conventions 

Our metric and spacetime index convention are de- 
fined by 77 M „ = diag(l, -1, -1, -1), x^ = (t,x,y,z) with 
fx = 0, 1, 2, 3. Whenever x appears alone, without indices 
attached, it always means the fi = 1 component of . 
The Einstein summation convention always applies un- 
less otherwise stated. The spacetime inner product of 
and 6 M is a- b; while a 2 = a^a^ . Bold fonts denote spatial 
vectors; for instance, a • b = iy-aV and a 2 = 5ija l a^ . 

Dirac action: The action in ([1]) define the dynamics of 
a Dirac fermion interacting with the photon vector po- 
tential A^. We make use of the Feynman slash notation. 
In particular, 

]pip = 7 /1 (9 (tt + ieA^ip 

Except in section (|V|, we employ the chiral basis for the 
{7^} matrices, 



which is, in turn, defined via the Pauli matrices, 



<J^ 




r" = (l 2 x2,^) 



(Al) 



1 

1 



-i 

1 



1 

-1 



Spinorial solutions: Except in section |Vj interaction 
picture perturbation theory will be used. It is carried 
out using the plane wave solutions to the non-interacting 
massive Dirac equation {i$ — m)^ = 0, where the posi- 
tive energy -0( fcs | + ) and negative energy ^(fc,^-) solutions 
with momentum k and spin ,s are 



V'(MI-) 



ik-. 



where 



(A2) 



V» = (12X2,-<T* 



The v 'a ■ k and \pd~k are the matrices a ■ p and a ■ p 
written in diagonalized form (i.e. UDU~~ l , where D is 
diagonal), with the eigenvalues replaced with their pos- 
itive square roots. By going to the rest frame of the 
particle, it can be seen that the and £ 2 are the spin 
up (s = 1) and spin down (s = 2) states for the fermion 
with respect to the basis of Pauli matrices used here. 
For the antifermion, they are the spin down (s' = 1) and 
spin up (s' = 2) states. These plane wave solutions are 
normalized such that 



4>\k s \±Ak',s'i±)d 3 x = (2tt) 3 S^ (k - k% s , (A3) 



States of definite hclicity are obtained by taking the 
spinors £ to be eigenspinors of a ■ k with eigenvalues 
±|k|. For the particle the positive eigenvalue corresponds 
to positive helicity and negative to negative; for the an- 
tiparticle the negative eigenvalue corresponds to positive 
helicity and vice-versa. 



